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Abstract 

We study the excitation spectra and the dynamical structure factor of quan- 
tum Hall states in a finite size system through exact diagonalization. Com- 
parison is made between the numerical results so obtained and the analytic 
results obtained from a modified RPA in the preceding companion paper [0. 
We find good agreement between the results at low energies. 
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I. INTRODUCTION 



By using a singular gauge transformation to convert electrons into fermions interacting 
through a Chern-Simons field, much progress has been made in understanding the physics 
of a two-dimensional interacting electron systems in strong magnetic fields at simple filling 
fractions [|T]-0]. At the mean field level, this transformation maps certain fractional quan- 
tized Hall states into integer quantized Hall states of transformed fermions so that some 
conventional many-body techniques may be used to tackle the problem. In the preceding 
companion paper an analytic procedure call the modified random phase approximation 
(modified RPA) has been constructed in the spherical geometry for evaluating the response 
of certain quantized Hall states to external perturbations. The purpose of this paper is 
to present numerical results of excitation spectra and the dynamical structure factors at 
selected filling fractions, calculated using exact diagonalization and the modified RPA, so 
that comparison can be made and the performance of the modified RPA can be evaluated. 
Specifically, we present numerical results from both exact diagonalization and modified RPA 
for systems of = 8 and A^ = 12 electrons at filling factor z/ = 1/3 and z/ = 3/7 respectively. 

As usual for numerical calculations, we work in the limit where the electron-electron 
interaction is weak compared to the cyclotron energy. For the exact diagonalization method, 
this means that all electrons are confined to the lowest Landau level. Within the modified 
RPA calculations, this means that we are working in the limit where the bare band mass 
rrib is infinitely small compared to the effective mass m*, whose scale is set by the electron- 
electron interaction strength. (In practice, we have carried out these calculations for the 
value m* /mi, = 50.) 

Our main purpose is to compare the exact and approximate calculations of the dynam- 
ical structure factor S{q, u) over a range of wavevectors, for the frequencies which are on 
the scale of the electron-electron energy, and therefore very much lower than the cyclotron 
frequency Uc- For such frequencies, the structure factor S{q, u) is equivalent to the projected 
structure factor S{q,u), in the exact calculation where all electron operators are projected 
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onto the lowest Landau level. To obtain additional information, however, we wish to com- 
pare various frequency moments of the projected structure factor with the corresponding 
quantities obtained from the modified RPA. (Specifically, we calculate contributions to the 
/-sum rule, the projected static structure factor S{q) , and the static wavevector-dependent 
compressibility.) For these comparisons, we define the projected structure factor in the mod- 
ified RPA by excluding the contribution of the "Kohn mode" , which occurs at the cyclotron 
frequency at g ~ 0. 

Our plan for this paper is as follows. In section ||, we describe the model and briefly 
discuss the method we use in the exact diagonalization. In addition, we review a few 



necessary points of the modified RPA calculation. In section |T|, we present results for the 
projected dynamical structure factor S{q, u) which describes the low energy response of the 
system to an external perturbation and gives the dispersion of the collective modes of the 
system. As well as showing S{q,u!), we will also calculate the dispersion of the collective 
mode in the single mode approximation (SMA), the contributions of the excitations in the 
lowest Landau level to the projected static structure factor S{q), the /-sum rule, and the 
compressibility sum rule. Our conclusions are summarized in section fV[ 



II. MODEL AND METHOD 

We will use the spherical geometry throughout this paper, so our system is given by a 
magnetic monopole of total flux A*",^ quanta at the center of the sphere with N electrons 
confined to the sphere's surface, interacting via a Coulomb potential. If the system were 
infinitely large (ie, a planar system), the filling fraction would be given by the ratio of the 
number of electrons to the number of flux quanta (z/ = -^). However, in a finite sized 
system, this relation is changed to 

N^=[u]-'N-X{iy). (1) 

where X(z/) is known as the "shift" of the state, and depends also on the topology of the 
system. In this work, we will be interested in filling fractions of the form Up = 2fp[- -^^^ 
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these states, in our spherical geometry, the shift is given by 



X{up) = 2 + p. 



(2) 



In the spherical geometry, by rotational invariance, the eigenstates of the system can 
be classified by the conserved quantum numbers (L, M), where L is the quantum number 
for the total angular momentum of an eigenstate and M is its z-component. In order to 
compare results on a finite sphere with those for an infinite planar system, we identify the 
wavevector q on the plane with the quantity L/R, where R is the radius of the sphere. This 
identification is exact in the limit R —>■ oo, but it is to some degree arbitrary for any finite 
system. (For example, the identification q — > \J L{L + 1)/R would be equally valid.) 

In keeping with the convention of the preceding companion paper we choose units 
of length such that the radius of the sphere is unity. The effective "size" of the sphere is 
then defined by the magnetic length Iq = S^^^"^ where S = is half the total number 

of flux quanta through the surface of the sphere. It is also convenient to choose the unit of 
energy to be e^/e/o which is the only energy scale in the problem when all of the electrons 
are confined to the lowest Landau level. 



Since exact diagonalization in the spherical geometry is a well established method, we 
will only describe it very briefly; more information can be found in Refs. and 0. The 
extraction of the projected dynamical structure factor S{q,uj) from such diagonalizations, 
however, is a method new to this work. 



are the monopole spherical harmonics Yf^^^fl), where I = S, S + 1, . . . , and m = —I, 
The corresponding eigenenergies are given by 



A. Exact Diagonalization 



As mentioned in the preceding companion paper [Q, the single particle energy eigenstates 




(3) 
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where ujc = hS/mf, is the cyclotron frequency, with mb the band mass of the electron. The 
lowest Landau level has a basis of the 25* + 1 states with / = 5*. With the single particle 
states given, ignoring the constant kinetic energy, the Hamiltonian of the many electron 
system, projected to the lowest Landau level, is written in second quantized notation as 

where Cj^ or C^rn creates or annihilates an electron in the single particle eigenstate Y;'^(r2), 
and Kn'j,m^;mi,m2 IS the matrix element of the spherical Coulomb interaction V{fli,Q,2) = 

where the vectors fl are restricted to the surface of the sphere. The relevant matrix elements 
are calculated more explicitly in Ref. 0. 

Now let us say a few words about the density operator in the spherical geometry. By 
definition, 

p(0) = ^5(f7-a), (6) 

i 

where fli is the position of the i*'* electron. Since p{fl) creates a neutral excitation when 
it acts on a state, it is useful to expand it in normal spherical harmonics, i.e. we define its 
Fourier components by 

piUn) = Y.nu^,), (7) 

i 

where Yi^m are the usual spherical harmonic (which are related to the monopole spherical 
harmonics p by Yim = ^m)- is easy to verify that pim has the expected properties: 

[L^,pim] = mpim, [L±,pim] = \J l{l + 1) - m{m ± 1) pim±i, (8) 

and if = 0, then 

i^Vim|$) = l{l + l)p«m|$), and L^pi„,\^) = mp;^|$). (9) 



In second quantized notation, we can write the density operator as 

p{n) =^^{n)^{n) (lo) 

where and are the operators that create and annihilate respectively an electron 

at the point ft. As usual, the eigenstate operators (Cj^^) can be easily related to the position 
space operators (ipK^)) 

cl^ = j dnYl{n)^\n), (ii) 

which can easily be inverted using the orthogonality property of the monopole spherical 
harmonics to yield 

i'\n) = Y.[Yi%mrci^. (12) 

l,m 

Thus we can write the density operator as 

pin) =Y. H Krn.m*yilrn,{^)ClrnA2,n... (13) 

li,mi l2,m2 

We can then take the angular momentum components to yield 

pl,m = J dnYi^^{n) p{fl,) = J2 H P(^"^,^l,"^l,^2,m2)C/^_„^C;2,m2 (14) 

where 

p(/,m,/i,mi,/2,m2) = / F,,™(^^)[F,f„^(^7)]*F,f,„^(^7) (15) 

which is evaluated explicitly in Appendix D of Ref It is then easy to project this operator 
to the lowest Landau level (Zi = ^2 = "S") to get 

Vim = Y.'Pih rn')Cl„^^^,Cs,m', (16) 

m' 

where 

/97TT / s s i \ ( s s i\ 

p{l,m,m') = i-ir^^-^-'{2S+l)f-^ . (17) 



An important property of at / = 1 is that 

Pi^o oc L^, pi 1 oc L+, and pi _i oc L_, (18) 

where the bar over the angular momentum operators indicate that they have been projected 
to the first Landau level also. This can be shown by writing the angular momentum operators 
in terms of eigenstate creation and annihilation operators p|, projecting to the first Landau 
level and comparing to the definition of the projected density operator. For example, we 
can write 

L, = j:^ClmCs,m (19) 
m 

which is easily shown to be proportional to 'Pi q once we have evaluated the 3-j coefficient in 
Eq. Using these properties (Eqs. [|18[) of p^ ^ we see that for |$o) in the first Landau 



level, L^|<I>o) = implies Pi,m|$o) = 0; ie the projected density operator cannot generate 
excited states at Z = 1 from a uniform state. 

Now we have come to the central quantity of our calculation. By definition, the projected 
dynamical structure factor is given by 

3(1, uj) = J2\{^^\piJ^o)\'Siiu + 60 - e,), (20) 

a 

which is independent of m by rotational invariance. By comparison, the full dynamical 
structure factor on a sphere is defined by 

S{1, u;) = E \{^a\pim\^o)\'S{u + eo - e^). (21) 

a 

In the limit where the band mass is taken to zero, while the frequency u and the electron- 
electron interaction are held fixed, the only energy states with finite — eg are states in 
which all electrons are restricted to the lowest Landau level. Thus if the frequency uj is fixed 
on the scale of the electron-electron interaction, the quantities 5'(/,ti;) and S{l,uj) become 
identical, as mentioned earlier. However, at energies on the scale of the cyclotron frequency, 
transitions between Landau levels will always be important, and the projected calculation 



will be inaccurate. In particular the projected structure factor will never display a Kohn 
mode and will never satisfy the /-sum rule. 

We should also mention that for the exact diagonalization it is possible to take advantage 
of the rotational symmetry of the spherical system to reduce the dimension of the relevant 
many-body Hilbert space. Specifically, we explicitly block-diagonalize the Hamiltonian by 
computing its matrix elements between many-body states with the same total angular mo- 
mentum and L^. To construct an eigenstate of the total angular momentum operator 
with eigenvalue /(/ + 1), we may in principle start from a randomly generated many-body 
state with the required and use the operator n^Y/l-^^ " ^'(^' + 1)] to project out the un- 
wanted components with other than l{l + 1). However, on a finite precision computer, a 
naive application of this method is highly unstable, i.e. the error introduced by the machine 
round-off grows exponentially. Fortunately, there are ways to stabilize the algorithm so that 
eigenstates of with up to machine precision can be obtained. Details will be published 
elsewhere |I0[]. 



B. Modified RPA 

As mentioned in the preceding companion paper [Q], there are in principle two free 
parameters in the modified RPA. They are the ratio of the effective mass to the band mass 

fe 2 2 

m*/mb, and a constant C = (;;^)/(^)) which is the ratio of the effective quasiparticle 
cyclotron energy to the Coulomb energy (which should be of order one). The estimate 
C ~ 0.3 has been made by Halperin, Lee, and Read by examining results from exact 
diagonalizations of small systems. We will use this value of C for all of our calculations in 
this paper. 

In order to project our results to the lowest Landau level, we should take the band mass 
to zero, or m*/mh oo. In practice, this is calculationally difficult, and we have actually 
used m*/m = 50 which is quite sufficient to separate the Kohn mode from the low energy 
excitations. We then use the prescription described the preceding companion paper |I| to 
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calculate the density-density response function Koo{L,u!) which is related to the dynamical 
structure factor via 

Smu;) = —lm[Koo{L,uj)]. (22) 

71 

As noted in [^1, the modified RPA is constructed to satisfy the /-sum rule, which one derives 
from the exact structure factor on the sphere. Within the units employed in this paper, this 
becomes 

r r o/r N L{L + l)Nm*C 
JO Stt nib S 

where the factor ^ (with S = ^ here) is simply a conversion of energy scales as described 
above. Note that this quantity diverges in the limit rrih/m* 0. The main contribution to 
this sum rule, however, is from the Kohn mode at (or above) the cyclotron frequency. All 
other modes described by the modified RPA are low energy excitations (on the energy scale 
set by the Coulomb interaction). Moreover, as discussed above, in the exact diagonalizations, 
the projection to the lowest Landau level restricts our results to only those modes on the 
energy scale of the Coulomb interaction. Thus, in order to "project" S{L,uj) calculated in 
the modified RPA to the first Landau level, we simply discard the high frequency Kohn 
mode and call the result the projected dynamical structure factor S{L,u). We then define 
a projected contribution to the /-sum rule by 

/(L) = / du uS{L,u). (24) 
Jo 

The value of f{L) remains finite for nib/m* — > for both the exact calculation and the 
modified RPA, so it is reasonable to try to compare the two. 

In a similar manner, we define the projected static structure factor S{L), and the con- 
tribution A(L) to the compressibility sum rule by 



S{L) = / dio S{L,uj), (25) 







AiL)=rdu^^. (26) 

Jo CO 
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In addition we shall report results for the average frequency entering in the single mode 
approximation (SMA) |p, which is simply defined by 

_ J,-'du;uJ{L,u;) _J{L) 
^ j;^dujS{L,u) S{L)' ^ ' 

The complete static structure factor S{L) is defined using S{L, uj) on the right hand side 
of Eq. (^) instead of S{L,uj). The contribution to S{L) from the Kohn mode, or from 
inter-Landau-level transitions does not diverge in the limit nih/m* 0, but it is a large 
finite contribution; therefore, to make sensible comparisons between the modified RPA and 
the exact calculations for the lowest Landau level, we use the projected structure factor in 
both cases. 

By contrast, in the limit mh/m* 0, the contribution to A(L) from the Kohn mode 
or from inter-Landau-level transitions is vanishingly small. Therefore, it would not matter 
whether we use S{L,u!) or S{L,u!) on the right hand side of Eq. (]26|). 



III. NUMERICAL RESULTS 

In this section we present and compare numerical results obtained from exact diagonaliza- 
tions and from the modified RPA, for two filling factors in the principal sequence Up = 2f^- 
The first set of data, shown in Figs. is obtained for a system of eight electrons at a 
filling factor z/ = 1/3 (As required by Eqs. ^ and 0, A'^ = 21 here). In the composite 
fermion picture, at the mean field level, the groundstate of the system corresponds to a 
completely filled lowest effective Landau level for the composite fermions. The second set 
of data, shown in Figs. |6|-p!0|, is obtained for a system of twelve electrons at a filling factor 



u = 3/7 {Nfj, = 23), where the groundstate of the system corresponds to the lowest three 
filled effective Landau levels of the composite fermions. 

As discussed in the Refs. and [^, in a planar system, we expect the modified 

RPA to be more accurate in the large p limit (ie, u approaching 1/2) where the motion 
of the composite fermions is semiclassical. In the present calculations however, we do not 
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necessarily expect the modified RPA to be more accurate for tlie = 3/7 = 3) state 
tlian for the = 1/3 {p = 1) state due to the problem of enhanced finite size effects for the 
u = 3/7 state. Although the u = 3/7 system contains more electrons (A^ = 12) than the 
u = 1/3 system (A^ = 8), it should be considered the "smaller" system in that the composite 
fermion magnetic length is much larger in the u = 3/7 system, and the number of electrons 
per Landau level is also smaller. Thus, any finite size effects will be accentuated in the 
z/ = 3/7 set of data, and we suspect that the modified RPA will be less accurate in this case. 

In Fig. ID, we show the projected spectral weights of S{L, u) for a system of eight electrons 
at a filling factor u = 1/3 {N^ = 21), calculated through exact diagonalization (upper graph) 
and using the modified RPA (lower graph). The weight of each delta function contribution 
to S{L, uj) at each (L, u) is proportional to the area of the corresponding shaded rectangle. 
Here and in the rest of this work, the energy scale is given in units of e^/(e/o)- As we 
discussed above, although there may be excitations at L = 1, these states all have zero 
weight in the projected dynamical structure factor. 

We see immediately from Fig. |I]that the projected dynamical structure factor calculated 
in the modified RPA agrees very well with the exact diagonalization at low energies. In 
particular, the distribution of spectral weights among the states in the lowest collective 
mode is very similar in the upper and lower graph. The small discrepancy in overall energy 
scales might be due to an inaccuracy in our guess of the energy scale conversion factor C 
described above. 

We note that the lowest collective mode dominates the spectral weight for angular mo- 
menta up to L = Lroton ~ 4 in the exact diagonalization, particularly near the roton 
minimum Lroton- In the modified RPA the lowest mode dominates only at smaller angular 
momenta. At higher angular momenta {L > Lroton), the exact diagonalization shows a fairly 
significant amount of spectral weight at energies E ~ 2Eroton- It is believed that the cor- 



responding states can roughly be represented as composite objects of two rotons [|TT|. We 
should not expect that such states would be properly modeled in the modified RPA since 
the RPA only represents single quasiparticle-quasihole excitations. However, as long as we 
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consider sufficiently low energy (not too much larger than E ~ 2Eroton), the modified RPA 
seems to agree very well with the exact diagonalization. 

At higher energies, in the continuum of the exact diagonalization spectra, we empirically 
observe that the spectral weights decay exponentially S{L,uj) ~ e'^^^ for a given L, where 
r ~ 0.03e^/e/o- In this high energy region, where the modified RPA differs significantly from 
the exact results, we see that the modified RPA yields a discreet spectrum of high energy 
modes with a large amount of spectral weight. The discreteness of the modified RPA is an 
obvious artifact arising from the neglect of processes where a higher energy mode can decay 
into several low energy modes with the same total energy and wavevector. The fact that the 
modified RPA seriously overestimates the total weight at high energies is a more significant 
limitation of the approximation. 

In Fig. we show unweighted excitation spectra for the same system of = 8 electrons 
at z/ = 1/3. The upper graph is the complete energy spectrum in the range < u < 0.45, 
< L < 11, from exact diagonalization regardless of whether or not the state contributes to 
S{q,uj). The lower graph is again the result of the modified RPA. Note that the dispersion 
relations and the overall energy scale of the lowest collective mode are very similar in the 
two cases. The black circles in the two graphs are the spectra of the single mode approxi- 
mation (SMA) computed from with Eq. [^, for the exact calculations and modified RPA 



respectively. We observe that in the exact diagonalization results, the SMA works well for 
angular momenta up to the roton minimum Lroton while in the modified RPA, the SMA fails 
at all but the lowest angular momentum. This failure of the SMA is due to the erroneous 
predictions of the modified RPA at high energy. In the modified RPA, as we increase the 
angular momentum, the main contribution to the spectral weight comes from a frequency 
that increases approximately proportionally to L{L + 1). This can be explained roughly by 
considering the modified RPA as a perturbation of a free electron system where it is well 
known that the excitations are confined to a band whose upper and lower energy boundaries 
both vary as g^. When a magnetic field is added to such a free electron system, a Kohn mode 
appears, and some of the low energy excitations can be changed, but we expect that most 
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of the high q excitations remains roughly in the same band. As discussed above, the single 
particle nature of the modified RPA allows these high energy modes to be given too much 
spectral weight. Thus we see how the the single mode approximation in the modified RPA 
is dominated by the poorly modeled higher frequencies in the modified RPA, thus giving 
very inaccurate results. 

In Fig. H, we show the projected static structure factor S{L) (see Eq. [^) for the 



same system of = 8 electrons at v = 1/3. The circles (solid line) are from the exact 
diagonalization whereas the triangles (dotted line) are the result of the modified RPA. The 
inset is the same modified RPA data on a reduced scale to more clearly show the behavior at 
higher angular momenta. Of course the projected static structure factor is only defined at 
integer values of L, and the lines (dotted and solid) are just guides for the eye. In the exact 
diagonalization (circles), we find that S{L) peaks around the roton minimum and decreases 
rapidly at higher angular momenta. We have studied its behavior and found empirically that 
it roughly obeys a Gaussian S{L) ~ Q~aL(L+i)^ where a is a constant. This reflects a similar 
empirically observed Gaussian decay of the matrix elements of the density operator at large 
angular momentum. The projected static structure factor in the modified RPA (triangles) is 
very similar to the exact result (circles) at small values of the angular momentum where the 
modified RPA does not predict extraneous high energy modes with large amounts of weight. 
At higher values of L the erroneous high energy modes cause a severe overestimation of the 
projected static structure factor. 

In Figs. ^ and |^, we show the contribution to the /-sum rule and compressibility sum 



rule as defined in Eqs. p4[ and p6[ for the same system of = 8 electrons at filling 
fraction z/ = 1/3. As above, the circles are from the exact diagonalization whereas the 
triangles are the result of the modified RPA, and the lines are guides for the eye. Once 
again the insets show the same modified RPA data on a reduced scale. In Fig. ^, we see 
that the contribution to the /-sum rule suffers from the same erroneous high energy modes. 
The compressibility sum rule (Fig. ^) is somewhat better estimated by the modified RPA 
because the contribution of the erroneous high energy modes are suppressed by a factor of 
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uo in this sum rule (see Eq. [0). This is somewhat encouraging, since the compressibihty 
sum is the most relevant of these quantities for describing the low energy behavior of the 
quantized Hall state. 

In Fig. 1^, we show the projected dynamical structure factor S{L,uj) for a system of 
= 12 electrons at a filling factor u = 3/7 {N^ = 23). Again, the spectral weights of 
S{L, uj) at each (L, u) are proportional to the area of the corresponding shaded rectangles. 
We observe again that the behaviors of S{L, u) from exact diagonalization (upper graph) and 
calculated in the modified RPA (lower graph) are semi-quantitatively similar. However, once 
again, the modified RPA results clearly overestimate the spectral weights at high energies. 

In Fig. ^ we show the corresponding complete excitation spectra of the system of 
A^ = 12 electrons at a filling factor u = 3/7. An interesting feature to note is that there 
are wiggles in the dispersion of the lowest collective mode in both the exact calculation 
and the modified RPA. Note particularly how the lowest energies at L = 3 and 5 lie above 
the lowest energies at L = 2,4, and 6. This behavior was observed previously in the exact 
calculation for this system size and filling factor by d'Ambrumenil and Morf [^]. These 
authors posed the question whether the even-odd alternation might be a spurious effect due 
to the finite size of the system. Based on the modified RPA analysis, however, we believe 
that the maxima and minima in the spectrum are genuine effects reflecting the maxima and 
minima that were previously found in the dispersion relation for the lowest excitation 
branch in the planar system using this approximation. More generally, for filling fractions 
of the form v = p/ (2p+ 1), at large values of p, it was predicted in Ref. that there should 
be a series of maxima and minima in the dispersion of the lowest excitation branch occuring 
respectively at wavevectors of the form q = 7r(n — \)/R* and q = T[{n + ^)/-R*, where n is a 
positive integer, and i?* is the "effective cyclotron radius" for the composite fermions. This 
quantity is given by the relation 

Rl = hk^/{eAB) = 2p/kp, (28) 

where AB is the deviation of the magnetic field from the field at z/ = 1/2, and the Fermi 
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wavevector kp is related to the electron density n^, by kp = (47me)^/^. Thus, at filling fraction 
I' = pI (2p +1), we have 

R: = y{2pf/\2p+lf'\ (29) 

If we now apply this formula to our system of particles on a sphere, with Iq^ = 
we find -R* ~ 2 at z/ = 3/7, with = 12. Thus one might expect maxima and minima 
to alternate with a period of AL = ir/R* ~ 1.6. Since the observed alternation with a 
period of AL = 2 is not far from this expectation, we believe it is appropriate to identify 
the oscillations as the expression in our finite size system of the oscillations predicted for 
the plane. In any case, the fact that this detail of the spectrum predicted from the modified 
RPA for the finite system is in good agreement with the exact calculations gives support to 
the belief that the oscillations predicted for the infinite system should also be present in an 
exact calculation. 

The black circles in Fig. ^ show the single mode approximation (Eq. [^) for u = 3/7. 
As we discussed in the case of u = 1/3, the single mode approximation clearly shows that 
the modified RPA predicts too much spectral weight at high energies. 

In Figs. H, ^ and |10|, we show the projected static structure factor (Eq. l2^), the 
contribution to the /-sum rule (Eq. |2^) and the contribution to the compressibility sum 
rule (Eq. [^) for the same system of A^ = 12 electrons at filling fraction u = 3/7. As above, 
the circles are always from the exact diagonalization and the triangles are from the modified 
RPA. Once again the insets show the modified RPA data on a reduced scale to show the 
behavior at higher angular momenta. Again, we see that the erroneous high energy modes 
in the modified RPA dominate the /-sum and the static structure factor at all but the lowest 
angular momenta. On the other hand, the compressibility sum (Fig. 0) is somewhat better 
represented in the modified RPA due to the suppression of these high energy modes. 
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IV. CONCLUSIONS 



We have studied the excitation spectra and the dynamical structure factor in quantum 
Hall states belonging to the principal sequence v = p/ (2p+ 1) through exact diagonalization 
and projected modified RPA. In particular, we have presented results from finite size exact 
diagonalization and from projected modified RPA on a system of = 8 electrons at = 1/3 
and a system of A^ = 12 electrons at z/ = 3/7. We find that the modified RPA works 
reasonably well at low energies, ie for energies ^ 0.2, in units of e^/e/o- (By comparison, the 
effective Fermi energy E^, = k^/{2m*) is approximately C/2 ^ 0.15, in these units.) 

By combining our results with previous analyses, we can arrive at some reasonable con- 
jectures for the excitation spectrum in a planar system, which should hold in particular for 
relatively large values of the parameter p. 

1. The present work supports the conjecture that the modified RPA gives a qualitatively 
correct description of the dispersion curve for the lowest branch of the excitation spectrum 
(quasiexciton mode) at u = p/ {2p + 1). The series of maxima and minima predicted by the 
modified RPA in the planar limit have analogs in the finite system which are also observed 
in the exact calculations. 

2. The general distribution of the spectral weight S{q, uj) for frequencies above the lowest 
excitation branch, up to energies of the order of ~ 0.2e^/ (e/o), for wavevectors in the range 
< < 2/cp is likely to be represented in a qualitatively correct fashion by the modified 
RPA. Of course certain details must clearly be wrong as has been previously discussed 0. 
Higher branches of the energy spectrum, which are undamped in the modified RPA, should 
actually be broadened into a continuum by decay into multiple excitations of lower lying 
modes. 

3. The modified RPA predicts that for large values of p, and wavevectors in the range 
kp/p <^ q ^ hp, the dominant contributions to the static compressibility sum rule should 
arise from frequencies u (x q^. We believe this prediction to be correct, but our finite size 
systems are too small to give any direct confirmation. 
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4. The modified RPA predicts that there should be sizable contributions to S{q,u;), for 
cu ^ 0.2e^/(e/o), for large wavevectors in the range where |g^ — 2m*ui\ ^ 2A;Fg. We believe 
these predictions to be spurious, however, as exact finite size system calculations show little 
weight in S{q,ijj) at high frequencies for any value of q. Of course the complete dynamical 
structure factor S{q, u) will contain additional contributions at very high frequencies arising 
from transitions between Landau levels. 

5. The above conjectures are of course predicated on the assumption that no instability 
or phase transition occurs at very large values of p, which would then invalidate our analysis. 
Finite size system calculations cannot rule out the possibility of such an instability at large 
P- 
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FIGURES 



FIG. 1. 

The weights of the projected dynamical structure factor S{L, cu) for = 8 electrons at 
filling fraction z/ = 1/3 (A'^^ = 21). The upper graph is obtained from exact diagonalization, 
and the lower graph is obtained from modified RPA. The magnitude of each delta function 
contribution to the projected dynamical structure factor S{L, u) is proportional to the area 
of the corresponding shaded rectangle. 



FIG. 2. 

The excitation spectra for N = 8 electrons at filling fraction u = 1/3. The upper graph 
is obtained from finite size exact diagonalization, and the lower graph is obtained from 
modified RPA. The black circles are spectra calculated in the single mode approximation 
(see Eq. 



FIG. 3. 

The projected static structure factor S{L) (See Eq. ||25|) for A^ = 8 electrons at filling 
fraction u = 1/3. The circles are obtained from finite size exact diagonalization, and the 
triangles are obtained from the modified RPA. The inset is the modified RPA data on an 
expanded scale to show the behavior at large L. The solid and dotted fines are guides for 
the eye. 



FIG. 4. 
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The contribution to the /-sum rule f{L) (See Eq. p4[) for = 8 electrons at filling 
fraction u = 1/3. The circles are obtained from finite size exact diagonalization, and the 
triangles are obtained from the modified RPA. The inset is the modified RPA data on an 
expanded scale to show the behavior at large L. The solid and dotted fines are guides for 
the eye. 



FIG. 5. 

The contribution to the compressibility sum rule A(L) (See Eq. |2^) for A^ = 8 electrons 
at filling fraction u = 1/3. The circles are obtained from finite size exact diagonalization, 
and the triangles are obtained from the modified RPA. The inset is the modified RPA data 
on an expanded scale to show the behavior at large L. The solid and dotted lines are guides 
for the eye. 



FIG. 6. 

The projected dynamical structure factor S{L, u) for A^ = 12 electrons at filling frac- 
tion z/ = 3/7 {N^ = 23). The upper graph is obtained from exact diagonalization, and the 
lower graph is obtained from modified RPA. The magnitude of each delta function contri- 
bution to the projected dynamical structure factor S{L, u) is proportional to the area of the 
corresponding shaded rectangle. 



FIG. 7. 
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The excitation spectra for = 12 electrons at filling fraction v = 3/7. The upper 
graph is obtained from finite size exact diagonalization, and the lower graph is obtained from 
modified RPA. The black circles are spectra calculated in the single mode approximation 
(see Eq. E^) 



FIG. 8. 

The projected static structure factor S{L) (See Eq. [^) for = 12 electrons at filling 
fraction u = 3/7. The circles are obtained from finite size exact diagonalization, and the 
triangles are obtained from the modified RPA. The inset is the modified RPA data on an 
expanded scale to show the behavior at large L. The solid and dotted lines are guides for 
the eye. 



FIG. 9. 

The contribution to the /-sum rule f{L) (See Eq. [Q) for A^ = 12 electrons at filling 
fraction u = 3/7. The circles are obtained from finite size exact diagonalization, and the 
triangles are obtained from the modified RPA. The inset is the modified RPA data on an 
expanded scale to show the behavior at large L. The solid and dotted lines are guides for 
the eye. 



FIG. 10. 
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The contribution to the compressibihty sum rule A(L) (See Eq. p6|) for = 12 electrons 
at filling fraction v = 3/7. The circles are obtained from finite size exact diagonalization, 
and the triangles are obtained from the modified RPA. The inset is the modified RPA data 
on an expanded scale to show the behavior at large L. The solid and dotted lines are guides 
for the eye. 
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